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Abstract 

We show that fluctuations of the fireball shape in the longitudinal direction generate nontrivial 
rapidity correlations that depend not only on the rapidity difference, y\ — y%, but also on the 
rapidity sum, y\ + 2/2 ■ This is explicitly demonstrated in a simple wounded nucleon model, and 
the general case is also discussed. We show how to extract different components of the fluctuating 
fireball shape from the measured two-particle rapidity correlation function. The experimental 
possibility of studying the longitudinal initial conditions in heavy-ion and proton-proton collisions 
is emphasized. 
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I. INTRODUCTION 



When a fireball produced in heavy-ion collisions is studied in an azimuthal angle, it is 
very useful to parametrize its initial shape with the help of the Fourier decomposition [1-3]. 
This decomposition clarifies trends in the elliptic flow, the fluctuations in elliptic flow, and 
the higher harmonics flows, which are associated with the triangularity and other shape 
parameters of the fireball [4, 5]. 

The harmonic analysis gives important information about the transverse initial conditions 
in heavy-ion collisions, and the mechanism of the subsequent evolution of the produced 
fireball - see e.g. [6]. In particular, the success of the hydrodynamic model in describing the 
v n data at RHIC and LHC places new constraints on the initial conditions in the transverse 
direction. 

A similar idea can be applied to study the shape of the fireball in the longitudinal direc- 
tion. In this paper we will focus on rapidity, y, but our arguments hold for any longitudinal 
variable. This study was initiated in Ref. [7] (see also [8]) where it was argued that long- 
range rapidity correlations can be interpreted in terms of the fluctuating rapidity density of 
the created fireball. When applied to the STAR data [9], a significant asymmetric component 
in the fireball's rapidity shape was found in the most central Au+Au collisions. 

In this paper we extend the discussion presented in Ref. [7]. We demonstrate that the 
fluctuations in the fireball rapidity density result in a nontrivial structure of the rapidity 
correlation function, and propose to study the additional components beyond asymmetry 
described above. The experimental method to extract various components is also discussed. 

The structure of this paper is following. In the next section we discuss the problem in a 
simple model. We show that an event-by-event difference between the number of wounded 
nucleons in the target and the projectile results in a long-range asymmetry of the fireball. 
We derive the correlation function and show that it depends on both the rapidity difference 
and the rapidity sum. In section 3, we decompose the different components of the fireball 
rapidity density with Chebyshev polynomials, and show how to extract the strength of these 
components from the measured rapidity correlation function. The practical application of 
this idea is discussed in section 4, where we also include several comments. We summarize 
our paper with the conclusions in section 5. 

II. SIMPLE MODEL 

In this section we explicitly demonstrate in a simple model that an event-by-event global 
fluctuations of the fireball rapidity shape lead to a nontrivial two-particle rapidity correlation 
function. 

For a given heavy ion event, we denote the number of wounded nucleons moving to the 
left and to the right with wl and wr respectively. The average over many events will be 
denoted by (wl) and (wr). In collisions characterized by (wl) ^ (wr), the single particle 
rapidity distribution is asymmetric with respect to y — 0, where y represents the rapidity in 
the center-of-mass frame. This asymmetry is clearly evident in d+Au collisions as measured 
at RHIC [10], and is easily reproduced by practically all models of heavy- ion collisions - see 
e.g. [11-14]. 
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In symmetric heavy ion collisions (Au+Au for example) we have (wl) = (%), and 
the single particle rapidity distribution is obviously symmetric with respect to y = [15]. 
However, this distribution is symmetric only when averaged over many events. In a single 
event the shape (in rapidity) of the fireball may be asymmetric since wl 7^ wr. 1 Indeed, in a 
single event the number of wounded nucleons going to the left may differ from the number of 
wounded nucleons going to the right. As discussed below, the asymmetry can be quantified 
by ((wl — %) 2 }, which is significantly larger than zero in Au+Au collisions. 

It is a useful exercise to calculate in a simple model the two-particle rapidity correlation 
function originating from fluctuations in wl — Wr- Here we consider the wounded nucleon 
model [16], which is a very useful model for understanding many features of heavy- ion 
data [11, 17, 18]. To simplify our considerations, let us assume that the single particle 
rapidity distribution measured in d+Au collisions can be approximated by a liner function 
of rapidity 2 . Consequently the distribution from a single wounded nucleon is also a linear 
function of rapidity. In the wounded nucleon model, the single particle distribution at a 
given wl and wr is given by [17] 

p(y;w L ,w R ) = w R (a + by) + w L (a-by) 

= a(w L + w R ) - by (w L - w R ) , (2.1) 

where a + by is the rapidity distribution from a right-mowing wounded nucleon, and a — by 
is the contribution form a left-mover. As seen from above equation we have an asymmetric 
component that is proportional to y. Assuming further that at a given wl and wr there are 
no correlations in the system 3 , the two-particle rapidity distribution at a given wl and wr 
is 

P2(Vi,V2;wl,wr) = p(yi;w L ,WR)p(y 2 ;w L ,w R ) 

= a 2 (w L + w R f - ab(w 2 L - w 2 R )(y! + y 2 ) + yiy 2 b 2 (w L - w R ) 2 . (2.2) 

Summing Eq. (2.2) over wl and Wr with an appropriate probability distribution, P(wl, wr), 
we obtain the experimentally accessible two-particle rapidity distribution. Taking (w L ) = 
(w R ), corresponding to symmetric Au+Au collisions, we obtain 

Pi{y\,y-i) = a 2 ((w l + wr) 2 ) +yiy 2 b 2 {{w L -w R ) 2 ) . (2.3) 

Consequently, the two-particle rapidity correlation function reads 

C(yi,y 2 ) = 0j(j/i,J/2) ~ p(yi)p(V2) , 

= a 2 [(wl) - (w + ) 2 } + y lV2 b 2 (w 2 _) , (2.4) 

where w + = wl + wr and W- = wl — wr. As seen from Eq. (2.4) the fluctuations in wl — wr 
result in a nontrivial rapidity structure of the two-particle correlation function. C(yi,y 2 ) 



1 The asymmetry due to the finite number of produced particles is not relevant to this analysis. 

2 This assumption is quite reasonable outside the fragmentation regions [10]. 

3 We want to study correlations originating only from shape fluctuations and we neglect all other correla- 
tions. We will come back to this point in section 4. 
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FIG. 1: Asymmetry in wounded nucleons in p+Au (left) and Au+Au (right) collisions at \fs = 200 
GeV as a function of the impact parameter b. 

depends not only on the rapidity difference, y- = y\ — y 2 , but also on the rapidity sum, 
V+ = 2/1 + 2/2- Indeed, the correlation function 



decreases as a function of rapidity difference, y_, and increases as a function of rapidity 
sum, y + . This dependence on y + can distinguish fluctuations of the fireball shape from well 
known sources of correlations (such as resonance decays) that depend mainly on y\ — y 2 . 

Eq. (2.4) should be taken as an illustration of the problem we would like to adress in 
this paper. Despite its simplicity, the model result shows quite convincingly that the event- 
by-event asymmetry of the fireball shape that is present in symmetric heavy-ion collisions 
can lead to interesting rapidity correlations [7]. Obviously there can be more complicated 
sources of this asymmetry in more realistic models, e.g. the difference in the number of 
flux-tubes in the CGC/Glasma approach [19, 20]. 

In Fig. 1 we present ((wl — w e) 2 ) divided by the total number of wounded nucleons 
(wl + Wr). We performed our calculations at y/s = 200 GeV in Au+Au and p+Au col- 
lisions in the Monte-Carlo Glauber model. It is interesting to notice that ({wl — wr) 2 ) 
in Au+Au collisions is quite large and comparable in magnitude to the total number of 
wounded nucleons. Thus, event-by-event rapidity fluctuations are of order ~ 1/ (wl + wr) 1 / 2 
compared to the average, and this is large compared to normal statistical fluctuations of 
order ~ 1/7V 1 / 2 , where N is the number of particles. 

In the next section we generalize Eq. (2.4) to analyze arbitrary rapidity fluctuations of 
the fireball shape. 

III. GENERAL SHAPE FLUCTUATIONS 

In the previous section we discussed the asymmetric component of the fireball rapidity 
density, originating from a non-zero values of wl — Wr- As seen in Eq. (2.1) for the simple 



C{yi,y 2 ) ~ 2/12/2& 2 (wt) 



^ 2 (2/i-2/^)<^>, 



(2.5) 
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model of the previous section, the single particle rapidity distribution at a given wl — wr 
is proportional to rapidity y, i.e. the fireball is denser on one side of the rapidity window 
than on the other. There may be different sources of this asymmetry such as the left-right 
difference in the number of collisions, or the difference in the number of asymmetric long- 
range flux-tubes in the CGC/Glasma approach [19, 20]. Let us denote the parameter that 
characterizes this asymmetry by a\. In our simple model, a% oc Wi — Wr. 

Equation (2.1) also contains a term that is proportional to the total number of wounded 
nucleons. Fluctuations of this quantity lead to symmetric, rapidity independent, fluctuations 
of the whole fireball. This can naturally originate from impact parameter fluctuations, which 
are always present in heavy-ion collisions. Let ao denote the parameter that characterizes 
this source of fluctuation. In our simple model, a oc wl + wr. 

The natural question arises if there are more components in the fluctuating shape of the 
fireball. For example, a "butterfly" component would characterize a symmetric fireball with 
higer (or lower) density on both sides of the midrapidity region, and lower (or higher) density 
at mid-rapidity. 4 The single particle rapidity distribution affected by this component would 
be proportional to y 2 . Let us denote by a<i the parameter that characterizes the strength of 
this effect. We will parametrize the fireball asymmetry and the butterfly component with 
the two Chebyshev polynomials, T\{y/Y) and T2(y/Y), which are shown in Fig. 2. 

Similarly, we can introduce additional components to fully parametrize shape fluctuations 
in rapidity. Thus, it is tempting to expand the single particle rapidity distribution at a given 




FIG. 2: Two components of the fluctuating fireball rapidity density: An asymmetry (solid black 
line) and a butterfly (dashed blue line). 



4 Such "butterfly" fluctuations are suggested by the measured forward-backward rapidity correlations at 
RHIC [9]. For a fixed number of particles at midrapidity, it was observed that the particle yields in pair 
of narrow rapidity bins located symmetrically about midrapidity strongly fluctuate. Surprisingly, these 
fluctuations are also strongly correlated. Thus, even if the density is approximately fixed in the middle of 
the fireball, both sides of the fireball fluctuate together. The physical origin of this correlation is currently 
under investigation - see Refs. [21, 22]. 
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ao, ai, ... in terms of the orthogonal polynomials 



p(y; do, ai,...) = p(y) 1 + V" . a(Ti (y/y) , (3.1) 

'1=0 



where p(?/) is the single particle distribution averaged over a , a±, Here we have expanded 

the distribution in Chebyshev polynomials 5 , but other choices are certainly possible. The 
parameter Y characterizes the scale of long-range rapidity fluctuations in the system. We 
will discuss reasonable values of Y in the next section. Averaging both sides of Eq. (3.1) 
over a ,ai, ... with an appropriate probability distribution, P(a ,ai, ...), we obtain (a*) = 
for all i > 0. 

Assuming that at a given ao,a±, ... there are no other large sources of long-range rapidity 
correlations, the two-particle rapidity distribution is 

Pi(yi, V2\ ao, ai, •••) = p{yi\ a , oi, --)p{y2] a , a u ...) . (3.2) 

Taking an average over a; and subtracting p(yi)p(y2), we obtain the two-particle rapidity 
correlation function 

C( yi ,y 2 ) = p( yi )p(y 2 ) [V (a l a k )T l (yi/Y)T k (y2/Y)] . (3.3) 

It is useful to recall the physical meaning of the first few terms in Eq. (3.3): (oq) 
represents the rapidity independent fluctuations of the fireball as a whole, (aocti) y 2 de- 
scribes the correlation between rapidity independent fluctuations of the fireball and its 
asymmetry, (af) y\y 2 is the asymmetric component discussed in the previous section, and 
(a|) [2 (|/i/y) 2 — 1][2 (2/2/y) 2 — 1] represents the butterfly contribution, etc. . 

From the previous section, we know that the asymmetric component, (af), introduces a 
long-range rapidity correlation that is a decreasing function of the rapidity difference, y_ = 
Hi ~ H2-I and an increasing function of the rapidity sum, y + = y\ + y 2 . It is a straightforward 
to verify that the rapidity structure originating from the butterfly component leads to a 
correlation function that is decreasing both in y_ and y + . 

To conclude this section, we point out that the values of (a^afc) can be extracted directly 
from the correlation function C(yi,y2)- Using the orthogonality of Chebyshev polynomials 

Ti(x)T k (x) (1 - x 2 )~ 1/2 dx = a5 i>k , (3.4) 

where Co = vr and q = ir/2 for i > 0, we obtain 

1 f Y C( Vl ,y 2 ) T t ( yi /Y)T k (y 2 /Y) d Vl dy 2 

{aiClk) c l c k J. Y p(y 1 )p(y2)[l-(y 1 /Yyf 2 [l-(y 2 /Yyf 2 Y 2 ' K ] 

In the next section we will discuss how Eq. (3.5) can be used in practice. 



For example: T (x) = 1, T x {x) = x, T 2 (x) = 2x 2 - 1, T 3 (x) = 4t 3 - 3a; etc. 
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IV. COMMENTS 



In this section we list several comments to clarify the analysis presented in this paper. 
In deriving Eq. (3.3) and Eq. (2.4), we assumed that at a given ao,a\,... there are no 
correlations in the system. In other words, the only sources of correlations are fluctuations 
in the fireball rapidity density. Unfortunately, short-range correlations may contribute to 
the left-hand side of Eq. (3.3), and contaminate the signal coming from the shape fluctu- 
ations. Particularly problematic may be the correlations from resonance decays and local 
local charge conservation [23, 24]. These problems can be mitigated by studying Eq. (3.3) 
for positive and negative particles separately, which significantly reduces these unwanted 
backrounds. Moreover, the dependence of the correlation function on the rapidity sum, 
y + , can be used to distinguish between rapidity density fluctuations and the short-range 
correlations of the background. 

One could also worry that at a given Wl + wr the distribution of final particles is given 
approximately by a negative binomial distribution (NBD) [25, 26] which introduces long- 
range rapidity correlations into the system. This concern is unfounded, however, because 
the NBD leads to the following two-particle rapidity distribution 6 



where k measures deviation from Poisson distribution. As seen from Eq. (3.3), the NBD p 2 
influences only (oq). In fact, this is the expected result, since the NBD can be viewed as a 
rapidity independent fluctuation of the whole fireball. 

As pointed out in the previous section, it is not obvious what is the appropriate value of Y 
in the preceding formulas. Clearly, Y parametrizes the range of global rapidity fluctuations 
in the fireball density. For instance, at \fs = 200 GeV the single particle distribution in 
d+Au collisions [10] is approximately linear as a function of y for \y\ < 2. Thus, for the 
asymmetric component parametrized by a±, Y « 2 is a reasonable choice. For higher and 
lower energies, this parameter can be rescaled by the ratio of beam rapidities. This value 
of Y roughly corresponds to the size of the thermal fireball, and it is plausible that higher 
components, if they exist, are present in this region. If the measurement is performed in the 
smaller window than [— Y, Y], fitting the measured correlation function with Eq. (3.3) can 
determine the parameters of interest, (ajOfc). 

Another choice is to assign Y to the rapidity interval of the measurement, and to investi- 
gate how the components (ajOfc) change when this rapidity scale is varied. It is possible that 
the different fluctuating components are visible at different rapidity scales, and a systematic 
study of this sort can sort out these differences. 

It would be very interesting to compare the strengths of different components, (ajOfc), 
between heavy-ion and proton-proton collisions. This could reveal interesting differences in 
the longitudinal initial conditions between these two systems. For instance, in Ref. [7] it 
was shown that the asymmetric component is significantly stronger in p+p collisions than 
in central Au+Au collisions. 

Further, the ideas presented in this paper could be extended by incorporating the multi- 
bin analysis proposed in Refs. [27-29]. This multi-bin analysis can be used to investigate 



6 We sample particles from NBD and distribute them randomly in rapidity according to p(y). 




(4.1) 
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the different sources of particles production, providing a detailed picuture of the fireball in 
the longitudinal direction. Finally, we point out that the results obtained in this paper can 
be easily generalized to three- and many-particle correlation functions. 

V. CONCLUSIONS 

In conclusion, we showed that event-by-event fluctuations of the fireball rapidity density 
introduce interesting rapidity correlations that depend both on the rapidity difference, y\ — 
?/2, and the rapidity sum, y\ + y%. We demonstrated this explicitly in the wounded nucleon 
model, where an event-by-event difference between the number of wounded nucleons in a 
target and a projectile, wi — wr, leads to the long-range asymmetry of the fireball. The 
resulting correlation function in symmetric A+A collisions is given in Eq. (2.4). 

We further proposed to expand the measured two-particle rapidity correlation function in 
a series of the Chebyshev polynomials (see Eq. (3.3)), where each polynomial represents a 
different component of the fireball's fluctuating rapidity density. The quadratic polynomial 
in this expansion describes the "butterfly" fluctuations described above, which are suggested 
by recent measurements at RHIC. The coefficients of this expansion, (aj%), characterize 
the strength of various components, and we propose to extract these coefficients from the 
measured correlation function. This can reveal nontrivial information about the structure of 
the fireball in the longitudinal direction, and can test various models of particle production 
in hadronic collisions. 
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